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, Abstract. Quasi-local variables, i.e. quantities whose values can be derived from 

physics accessible within an arbitrarily small neighborhood of a spacetime point, are 
used to construct the equation of state (EoS) for the anisotropic fluid in spherical 
O ! symmetry. One parameter families of equilibrium solutions are obtained making it 

possible to assess stability properties by means of the standard M(R) method. Normal 
£h . modes of radial pulsation are computed as well and are found to confirm the onset of 

instability as predicted by the M(R) method. As an example, a stable configuration 
with outwardly increasing energy density in the core is obtained with a simple quasi- 
local extension of the polytropic EoS. It is also found that the loss of stability occurs 
QO ' at higher surface compactness when the anisotropy of pressures is present 
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PACS numbers: 04.40.-b,04.40.Dg 
1. Introduction 

The most general spherically symmetric solution to Einstein equations allows for the 
anisotropy of principal pressures; the radial pressure may differ from the transverse 
pressure. While the isotropic fluid with specified equation of state (EoS) is the most 
common source for modeling Newtonian or relativistic stars, anisotropy of principal 
pressures can be found when the source is derived from field theories, e.g. when scalar 
fields are considered as in boson stars pQ, or in configurations involving electrically 
charged matter. One also finds that anisotropy of pressures is present in the so called 
"exotic solutions" to Einstein equations such as wormholes [21 [3] or gravastars [U [5]. 
Interesting features of relativistic anisotropic spheres were pointed out as early as 1933. 
by Lemaitre [6j [7] , while ref . [8] is widely considered as the beginning of the epoch of 
more active research, see [9] for a review. Recently, general algorithms for generating 
static anisotropic solutions to Einstein equations have been formulated (see [lOj [TT] 
and references therein.) However, whether such solutions are considered more or less 
physically acceptable depends on a number of criteria which include fulfillment of various 
energy conditions of general relativity, the requirement that the speed of sound in the 
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fluid is subluminal, and perhaps most importantly, stability of the configurations with 
respect to perturbations. 

The most notable technique for the analysis of stability of isotropic relativistic stars 
with respect to radial perturbations is that of Chandrasekhar [12J. In this procedure 
the time dependent Einstein equations are linearized around the equilibrium solution to 
yield a linear wave equation for radial perturbations, which together with appropriate 
boundary conditions comprises an eigenvalue problem of the Sturm-Liouville type. 
Positive eigenvalues are interpreted as (squares of) frequencies of normal modes of radial 
pulsations, while negative eigenvalues imply exponentially growing perturbations, i.e. 
instability of the sphere (for in-depth coverage see [T3l [T^] ). 

Another field-proven technique, known as the M(R) method, can be applied to 
locate the onset of instability in one-parameter families of equilibrium solutions to 
Einstein equations satisfying the same EoS of cold matter [131 US]- Using the central 
energy density as the parameter and increasing its value, the first maximum in the 
M(R) curve indicates the boundary between the stable and unstable configurations. 
While more limited both in its scope and insights into the dynamics of the system it can 
provide relative to the full analysis of the normal modes, the M(R) technique is very 
simple in circumstances where for the given EoS the equilibrium configurations can be 
generated over a sufficiently wide parameter range. 

In the context of anisotropic stars the normal modes technique was applied to study 
the stability of some specific configurations in [16J and also [171 [18]. However, the M(R) 
method could not be straightforwardly applied because most of the solution-generating 
techniques mentioned above could not generate sequences of equilibrium configurations 
corresponding to one fixed EoS. The only framework outside of isotropic spheres where 
the M(R) method was successfully applied is that of elastic stars. In the second paper 
in the series beginning with [T9l [20] the M(R) method was shown to predict instability 
exactly at the configuration that has vanishing frequency of the fundamental normal 
mode of radial pulsation, i.e. the two methods were shown to give compatible results. 

In this paper we will first show how a sequence of anisotropic spheres of different 
masses and surface radii can be generated from one given EoS. This will enable us 
to use the M(R) method to study stability, and we will compare the results to those 
of the analysis of normal modes. The paper begins (section [2]) with the discussion 
of the concept of the quasi-local EoS for the anisotropic fluid in spherical symmetry. 
In section [3] we derive the linear wave equation for radial perturbations for a specific 
quasi- local EoS where the radial pressure is described by the arbitrary barotropic EoS, 
P = p{p)i while the anisotropy is taken to be the arbitrary function of the energy density 
and the compactness, \x = 2m/r, which is a quasi-local variable. In section [4] we apply 
the above concepts to specific examples; we study the structure and stability of solutions 
obtained from the quasi-local EoS consisting of the polytropic EoS governing the radial 
pressure and the anisotropy of pressures taken to be bilinear in the radial pressure and 
the compactness. We close with a brief discussion of the results in section El 
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2. The equation of state 

Energy-momentum tensor of the perfect fluid can be written asT = pu®u + p(g + 
u ® u), where p is the energy density and p is the isotropic pressure of the fluid, u is its 
four- velocity normalized so that u-u = — 1 and g is the metric tensor. The combination 
g + u <g) u is the projection tensor onto the spatial hypersurface orthogonal to the fluid 
four-velocity. The notion of the anisotropic fluid allows the pressures to differ among 
spatial directions. In particular, in spherical symmetry, the radial pressure which we 
denote by p may differ from the transverse pressure which we denote by q. To write 
down the energy-momentum tensor of the anisotropic fluid in spherical symmetry, in 
addition to the four-velocity u, one introduces the 'radial' spatial unit vector k that is 
orthogonal to the four-velocity and has no angular components, i.e. k • k = 1, k • u = 0. 
The energy-momentum tensor of the anisotropic fluid can now be written as 

T = pu<g)u + pk<g)k + g(g + u<g}u-k®k) (1) 

where g + u®u — k ® k is the projection tensor onto the 2-surface orthogonal to both 
u and k. Since the vector k is not defined at the centre of symmetry the anisotropy of 
pressures must vanish at this point. 

To illustrate the role of the anisotropy of pressures in the dynamics of the fluid it 
is worthwhile to consider the first law of thermodynamics: 

d(pV) = -dW + TdS, (2) 

where V is the co-moving 3- volume of the fluid element, p and T are respectively the 
energy density and temperature of the fluid, S is the entropy and W is the work done 
by the fluid element. Assuming no energy flow among the neighboring fluid elements we 
have dS = 0. In the case of the isotropic pressure the work done by the fluid element 
due to the change in its volume is dW = p dV. Dividing by the element of the proper 
time along the world line of the fluid element one can write 
dW dV Tr/ „ , 

- d y=P^=P V (V-u) (3) 

where the relative change of the volume, (dV/dr)/V = V • u, can be recognized as the 
spacetime expansion scalar. Applying the chain rule on the lhs of (j2J), dividing by dr 
and using ([31) , for the isotropic fluid one obtains the well known relation 

u ■ Vp = -(p + p) V • u. (4) 

Turning now to the work done by the element of the anisotropic fluid we must, 
due to the difference among the principal pressures, accordingly distinguish between 
different directions of expansion. In spherical symmetry we distinguish between the 
radial and the transverse expansion: 

V • u = (V • u) rad . + (V • u) tr ., (5) 

where the radial expansion can be written as (V • u) ra( j. = k • (k • Vu). The work due 
to the expansion of the element of the anisotropic fluid is 

dW = pdV Iad . + qdV tv . = (pV (V ■ u) rad . + qV (V • u) tr .)dr, (6) 
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and the anisotropic analogue of (jl]) follows as 

u • Vp = -p V • u - p ( V • u) rad . - q (V ■ u) tr . 
= ~(P + ?)V • u - (p - g)(V ■ u) rad . 
= ~{P + <?) v • u - (p - g) k • (k • Vu). 



(7) 



The last term on the rhs is proportional to the anisotropy of the pressures and vanishes 
in case of the isotropic fluid, thus reducing (JTJ) to (SJ). However, (I7j) can also be obtained 
from the energy-momentum tensor of the anisotropic fluid ([1]) by taking the component 
of the conservation law VT = directed along the world line of the fluid element, 
u VT = 0, while the projection of VT = onto the hypersurface orthogonal to the four- 
velocity gives the anisotropic version of the Euler equation of relativistic hydrodynamics. 

The equation of state (EoS) of the isotropic fluid is usually understood as a relation 
among the energy density and the pressure of the fluid, and possibly other local fluid 
variables such as temperature, baryon density, entropy per baryon, etc. These variables 
are considered local because their values refer to the state of the fluid at the particular 
spacetime point, rendering the EoS independent of any non-local information. The local 
EoS describing the isotropic fluid can be cast in the form of a single equation: 



dots representing additional local variables. The simplest case is the barotropic EoS 
f{PiV) — 0, or as it is often written p = p(p), that can be taken to, e.g., describe the 
stellar matter at the end of nuclear burning and at zero temperature. It is well known 
[2Tj that the barotropic EoS suffices to close the static Einstein equations in spherical 
symmetry, and as shown in [22], a given barotropic EoS with pressure monotonically 
increasing in energy density guarantees the existence of a family of solutions with central 
energy density as the sole parameter. 

In analogy to the case of the barotropic EoS for the isotropic fluid, one could expect 
the simplest anisotropic EoS to be of the form p = p(p) for the radial, and q = q(p) for 
the transverse pressure. However, a simple argument can be given to show that such 
EoS is overly restrictive to allow for a family of solutions with different central energy 
densities. We begin by assuming the existence of a static spherically symmetric solution 
governed by the EoS of the form p = p(p) and q = q(p), with energy density po at r = 
(centre of symmetry) and non- vanishing anisotropy of pressures at some r = r\ > 0. 
Since at r = the anisotropy of pressures must vanish we have p(po) = <?(po)- At r = r\ 
we have p(pi) 7^ l{pi) by assumption, so it follows that pi ^ pq. If we now consider the 
possibility of constructing the solution with central energy density equal to pi, we find 
that the condition of vanishing anisotropy of pressures at r = cannot be satisfied. The 
EoS for the anisotropic fluid must, therefore, have more complex structure than what 
was just assumed. 

The idea behind the quasi-local EoS, which we borrow from (U [23], is to allow for 
the dependence on the quantities that can be derived from the geometry at a given 
point of the spacetime. In principle an arbitrarily small neighborhood of a spacetime 




(8) 
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point suffices to measure the (orthonormal frame) components of the Riemann tensor. 
These components, in turn, allow one to construct quantities that can be considered as 
quasi-local variables. In spherical symmetry the general quasi-local EoS will have the 
form of a system 

h(p,P,Q, ...;//,...) =0, f 2 (p,p,q,...;fi,...) = 0, (9) 

where p,p,q, . . . denote the local variables, and p, . . . denote the quasi-local variables. 
The quasi-local variables of special interest could be the curvature radius r or the 
compactness 2m(r)/r; for orthonormal frame components of the Riemann tensor in 
spherical symmetry see e.g. [3]. 

In contrast to the local or quasi-local variables as defined above, non-local variables 
were defined as arbitrary functionals over the whole, or some finite segment of the 
t = const, slice of the spacetime [23J [25]. While in principle arbitrary functionals can 
appear in the static limit of the EoS, this concept might turn out to be problematic 
when dynamics is considered. Namely, it is not clear how a perturbation of the fluid 
in the neighborhood of some spacetime point would affect the value of an arbitrary 
functional, and consequently the state of the fluid, elsewhere in the spacetime. The 
'immediate effect' of the perturbation on the functional would violate the principle of 
causality, while physically acceptable (causal) propagation would require the input of 
additional physics, implying that a non-local EoS alone does not fully describe the fluid. 
However, the 'average energy density up to r\ 

pW = tV/ 4r /2 7rp(r')dr', (10) 
3?""% J 

which is used as an example of a non-local variable in [24j|25] can, in fact, be understood 
as a quasi-local variable since p = 3m/4:r 3 iT = 3p/8iir 2 , where p = 2m/r and r are 
acceptable quasi- local variables [H 123]. 

3. The wave equation for radial perturbations 

Here we derive the linear wave equation for the radial perturbations of the anisotropic 
fluid governed by the specific form of the quasi- local EoS We assume that in the 
rest frame of the fluid its energy density p and the pressures p and q are related by 

P = p(p), q-p = a(p; fi), (11) 

where the quasi-local variable p is the compactness, p = 2m/r, m being the usual mass 
function. The procedure that follows is a generalization of the one originally carried 
out for the isotropic case by Chandrasekhar [12J, while our notation is (to the extent 
possible) compatible with Misner, Thorne and Wheeler [13], Ch. 26. 

The line element of the spherically symmetric spacetime in the standard coordinates 
t ,r, i? and (p is 

ds 2 = -e 2 * dt 2 + e 2A dr 2 + r 2 dfi 2 , (12) 
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where $ = $(t,r), A = A(t,r), and dQ 2 = di9 2 + sin 2 $ dip 2 is the line element on the 
unit 2-sphere. The metric function A is related to the compactness p, and to the mass 
function m by 

fi = 2m/r = 1 - e- 2A . (13) 

Denoting with £ = r) the radial displacement of a fluid element from its 
equilibrium position r at time t, the nonzero components of the four-velocity normalized 
so that u^u^ = — 1 are 

u t = ( e 2* _ e 2A e 2 )-V2 ; u r = (e 2* _ g^e 2 )- 1 /^, (14) 

where £ = <9 t £ = u r /u*. The nonzero components of the radial unit vector k^ satisfying 
= 1 and k^ = are 

= e-* +A (e 2 * - e 2A £ 2 )-V 2 £, k r = e*- A (e 2 * - e 2A (t 1/2 . (15) 

According to ([1]) the energy-momentum tensor is 

TV = (p + g)uV + + (p - (16) 

where p = p(t, r) is the energy density, p = pit, r) is the radial pressure and q = q(t, r) 
is the transverse pressure. 

One now assumes the existence of a static or equilibrium (time independent) 
solution to the Einstein equations and proceeds to study the dynamics of small 
perturbations. At a fixed point in the coordinate grid a quantity / = f(t,r) is 
decomposed as 

f(t,r)=f (r) + 5f(t,r), (17) 

where fo denotes the time independent equilibrium value of / and 5f denotes the so 
called Eulerean perturbation (in equilibrium 5f = 0.) In this way we decompose $, 
A, p, p and q, and through the rest of this section we omit the explicit notation of 
the dependence of variables on t and r; zero-subscripted equilibrium variables are 
understood to depend only on r, while 5-prefixed Eulerean perturbations depend on 
t and r. 

The Einstein equation, = SttT^, for the metric (|12p and the energy- 

momentum tensor ( |T6l) can now be expanded in powers of the displacement £ and 
the Eulerean perturbations 5$, SA, Sp, 5p and Sq. The zero-order terms govern the 
equilibrium configurations: 

Svrpo = r- 2 (-e- 2Ao (l - 2rA' ) + 1), (18a) 
8np = r- 2 (e- 2Ao (l + 2r%) - 1), (185) 
8tt g = r-V 2Ao ((r$ - rA' )(l + r$ ) + r 2 $ '), (18c) 

while the first order terms involve the fluid displacement and the the Eulerean 
perturbations: 

8tt 5p = 2r- 2 e- 2Ao ((l - 2rA' )5A + r 8A'), (19a) 

8n5p = 2r- 2 e~ 2A °(-(l + 2r$ )5A + r 5$'), (196) 
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8tt 5q = -e- 2 *°5A + r-V 2Ao ((1 + 2r& - rA' )r 5$' + r 2 o"$" 

-2((r$o -rAo)(l + r$o) + r 2 $' ')o"A- (1 + r$ / )r«JA / ), (19c) 
(Po + PoK = -2r- 2 e" 2Ao r 5A (19d) 

(overdot denotes d t , prime denotes d r .) 

From the EoS ( TTTT) one may derive the relations among the perturbation of the 
variables in the rest frame of the fluid. These are called the Lagrangean perturbations 
and are denoted by the A-prefix. From (II lj) we obtain 
dp 



Ap = -f-Ap, Aq - Ap 
dp 



da . da . 



(20) 



A fluid rest frame (Lagrangean) perturbation of a quantity Af is, for small 
perturbations, related to the perturbation of the same quantity seen in the coordinate 
frame (Eulerean perturbation) 5f with 

Af{t, r) = f{t, r + r)) - / (r) ~ Sf(r) + ti(r)Z(t, r). (21) 

Applying the above rule we express the perturbations of the fluid variables in terms of 
Eulerean perturbations, obtaining 

dp 



8p + 



dp 



(Sp + p'ot), 



dp 
dp 



+ 



da 
dp 



(5p + p' O + 2e 



-2A 



da 
dji 



(SA + A'S- 



(22) 
(23) 



The quantities in brackets follow from the assumed EoS and are to be evaluated at the 
equilibrium, as indicated by the zero subscript. 

The linearized Einstein equations (|18aH18c]) and ( jl9aH19a| . together with the 
perturbation of the EoS ( 122]) and ( 123|) . constitute a system of equations that yields 
the linear equation of motion for the fluid. After a lengthy algebraic procedure 
aimed at eliminating all Eulerean perturbations from the system so that only the fluid 
displacement £ and its derivatives remain, followed by the variable transformationr: 

C(t,r) = r 2 e-^ r k(t,r), (24) 
we arrive to the linear wave equation for the radial perturbations: 

C(t, r) = C 2 (r) C"(t, r) + d(r) (% r) + C (r) C(f , r). (25) 
The coefficients CAr) can be written as 



C 2 e 



d e 



-2*o+2A 



-2*o+2A 



C e 



-2$o+2A 



dp 
dp". . 

- $' + f 2# + A' 
%o - Po) 

K/°o +Po) 
(1 + r$ ) 2 + e 2A ° 



(26 a) 



dp 
dp". 



+ 



dp 
dp 



0+ 


dp 




2 


9a 


[dp_ 


.) 


r 


_<9p_ 



(266) 
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+ 



r 2 {po + Po) 



2r$ + (3(r$ - 1) + rA' 



dp 
d~p 



-3 + r 



2$' 



da 
dp 



+ 



4e 



-2 An 



r(po + Po) 



A' 



dp 
dp. 





da 
dpi 



dp 
dp" 



(26 c) 



The quantities in brackets must be computed from the EoS and evaluated using the 
equilibrium configuration variables. Following the standard procedure one assumes 
C(t,r) = u{r) exp(iwt) and the wave equation (|25|) reduces to the ordinary differential 
equation 

- u 2 u = C 2 u" + Cm' + C u, (27) 

which together with the appropriate boundary conditions imposed on u(r) (to be 
discussed below) constitutes the boundary value problem with eigenvalue u 2 . (The 
function u{r) should not be confused with the four-velocity u^.) Solutions with uj 2 > 
are oscillatory in time and therefore represent radial pulsations. Solutions with u> 2 < 
may exponentially grow in time and therefore imply instability of the equilibrium 
configuration with respect to radial perturbations. 

The boundary condition that we impose on u(r) at r = follows from the 
observation that, due to symmetry, the fluid element at r = cannot move. Therefore 
we have £(£, r = 0) =0 which with (1211) gives ((t, r = 0) = 0, and finally u(r = 0) = 0. 
Additionally, the Lagrangean perturbation of the energy density which can be written 
as 

e*° (1 
A P=~ — 



(28) 



-(?o-Po)C + (Po+Pq)C 
r" \r / 

must remain finite at r = from which it follows that u oc r 3 as r — > 0. The situation 
is more delicate at the surface of the star which we assume to be at finite radius r = R. 
There we require that the fluid displacement £(t,r = R) is finite at all times, which 
implies u(r = R) finite. Additionally, we must ensure that at r = R the Lagrangean 
perturbation of the radial pressure, Ap = [dp/dp]oAp, vanishes at all times as required 
by the junction conditions of the spherically symmetric spacetime with the exterior 
Schwarzschild spacetime [3]. It is important to note that if C 2 — > as r — > R then 
(T27|) is singular at r = R. Further analysis of the behaviour of u as r — > R, therefore, 
depends on the properties of the equilibrium configuration, or more fundamentally, on 
the properties of the EoS (see e.g. the discussion for the isotropic case in [13].) 
By multiplying ( )27l) with the weight function 

" C x {r)-C' 2 {r) 



W{r) = exp , ^ 
it assumes the Sturm-Liouville form, 



dr 



(29) 



- u 2 Wu = (Pu'Y + Qu, (30) 

where P = WC 2 and Q = WC . If W > over [0,R] and P, P', Q and W 
are continuous we have a regular Sturm-Liouville problem granting access to the 
main results of the Sturm-Liouville theory: if the eigenvalues are ordered so that 



Anisotropic stars with quasi-local EoS 



9 



u 2 < u)\ < . . . < uf < . . . < oo, then uf corresponds to the eigenfunction with i 
nodes in (0, R). This implies that in order to test the stability of the configuration with 
respect to radial perturbations it is sufficient to compute the eigenvalue u 2 corresponding 
to the solution without nodes. If Uq > then all solutions are oscillatory and the 
whole equilibrium configuration is considered stable, while if u 2 < it is considered 
unstable. (Solutions to the Sturm-Liouville eigenvalue problem are usually referred 
to as 'normal modes' due to the orthogonality of the eigenfunctions corresponding to 
different eigenvalues with respect to the weight function W(r).) 

We end this section by noting that if one takes a(p; p) = all the relations derived 
above reduce to the well known relations for the isotropic case (see e.g. [33] or [H], Ch. 
26.) We also note that instead of introducing the adiabatic constant T of the fluid as it is 
done in the standard procedure in the isotropic case where T = (p+p)p~ 1 dp/ dp we have, 
for clarity, preferred to work with the generic expressions throughout the procedure. 

4. Example: anisotropic polytropes 

To give an example where the concepts discussed so far are put to work in this section 
we study the properties of spherically symmetric solutions supported by the anisotropic 
fluid described with the quasi- local EoS: 

V = P(p) = kp 1+1/n , q-p = a(p; p) = ap(p) p. (31) 

The first relation is the polytropic EoS that is assumed to relate the radial pressure to 
the energy density. The parameter n > is known as the polytropic index. We have 
chosen this EoS because the solutions where the polytropic EoS governs the isotropic 
fluid, known as polytropes, are probably the most thoroughly studied self-gravitating 
objects in the literature, both in Newtonian and in the relativistic regime [261 [27]. The 
second relation in (13T1) is the ansatz for the anisotropy of pressures which we assume to 
be proportional to the radial pressure and to the compactness. The parameter a has 
the role of the anisotropy strength parameter. An ansatz similar to this one has been 
used to generate the gravastar solutions in [5]. An important feature of this ansatz is 
that the compactness which goes as r 2 when r — > ensures the required vanishing of the 
anisotropy of pressures at the centre of symmetry. Further motivation for this ansatz 
was drawn from astrophysical considerations that in the non-relativistic regime where 
the compactness is much smaller than unity the anisotropy of pressure is not expected 
to play an important role. Also, when combined with the polytropic EoS assumed for 
the radial pressure, this anisotropy ansatz ensures vanishing of the tangential pressure 
at low energy densities, e.g. at or near the surface of the star. 

The first part of our procedure is concerned with the equilibrium configurations. 
We begin by fixing the EoS by choosing a triple of parameters k, n and a in ( 13TI) . and 
proceed to compute a sequence of static (equilibrium) solutions to Einstein equations 
(H8aH18cj) corresponding to a sequence of central energy densities po(0) from a certain 
range. The equations to be solved for each po(0) can be compactly arranged as the 
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system of two coupled equations: 

m = 4vr r 2 p , (32) 

,2 (p +Po)Oo + 4vr r 3 p ) /qq\ 
Po = -{Qo-Po) z 77 , (33) 

where fl32|) involves the mass function mo which is related to the metric function A and 
compactness po by (fT3"j) , and fl3"3"|) is the anisotropic variant of the well-known Tolman- 
Oppenheimer-Volkov (TOV) equation (the zero-subscript indicates that the variables 
correspond to the equilibrium configuration and prime denotes d r .) By using the EoS 
to express po in terms of po, and by writing p' = [dp/dp] p' Q , TOV is converted into 
the differential equation for equilibrium energy density po(r). The solution to f[3"2"l) and 
(|33|) for the chosen EoS and central energy density po(0) provides us with the surface 
radius R and the total mass M = m(R) of the equilibrium configuration. Also, the r- 
dependence of m and p is known, from which p , q , p or A can be readily obtained. 
Finally, the metric function $o is obtained by integrating 

$' = m ° + 47Tr3p ° (34) 
r 2 (l-2m /r)' 1 ' 

which completes the calculation of the equilibrium solution! 

A sequence of equilibrium solutions corresponding to a given EoS, i.e. fixed triple 
n, k, a in (13T|) . can be generated by varying the central energy density Po(0) as the 
only parameter. Within the sequence the relation between the total mass M and the 
surface radius R can be studied. It can be represented as a curve in the M vs. R plane 
parametrized by the central energy density. According to the M(R) method [T5l fl3] . 
starting from the po(0) — > limit, all solutions with central energy density smaller than 
the one corresponding to the first maximum of the M(R) curve are stable with respect to 
radial perturbations. In general, each extremum of the M(R) curve represents a critical 
configuration in which some mode of radial pulsation changes its stability properties. 
Following the curve in the direction of increasing central energy density, at a critical 
point (extremum) through which the curve is making an anti-clockwise turn a mode 
loses the stability; in a clockwise turn through a critical point a previously unstable 
mode becomes stable. 

Our first example uses the EoS fixed by the triples n — 1, k — 100 km 2 , 
a = 0,±1,±2. This particular set of parameters is chosen so that for the isotropic 
case (a = 0) the M(R) curve reproduces the reference configuration of a relativistic 
polytrope presented in the table A. 18 of Kokkotas and Ruoff [28J. The M(R) curve 
corresponding to a = is shown with the thick line in the upper plot of figure [TJ The 
M(R) curves corresponding to positive and negative values of a show similar behaviour, 
i.e. they have a maximum which we expect to indicate the boundary between the stable 
and the unstable configurations. The second example shown in the lower plot of figured] 
uses n = 2, k = 5 km and a = 0, ±1, ±2. In both examples the maxima of the M(R) 
curves occur, as compared to the isotropic configuration, at higher surface compactness 
for a > 0. This hints that anisotropy of pressures with transverse pressure exceeding 
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Figure 1. Mass-radius curves for anisotropic stars with (radial) polytropic index 
n = 1 and k = 100km 2 (upper plot), and n = 2, k = 5km (lower plot): the M(R) 
curve for the a = (isotropic) configuration is shown with thick line, for a = 1,2 
with long-dashed, and for a = —1,-2 with short-dashed lines. Solid curves connect 
the configurations sharing the same central energy density and are drawn with step 
corresponding to factor 2 (central energy increases as one moves from right to left.) 
Circles indicate the maxima of the M(R) curves which coincide with w 2 — as 
obtained through the analysis of the normal modes. Straight thin lines originating 
from M = = R indicate the locus of points corresponding to surface compactness 
/ u s = 2M/ J R = l,2/3, 1 /2,2/5. 



the radial pressure (our a > 0) may be more efficient in supporting highly compact 
bodies against gravitational collapse than in the opposite case (our a < 0.) We can 
also observe that in our first example (n = 1) the maxima of the M(R) curves occur at 
lower central energy densities (as compared to the isotropic configurations) for positive 
a, while in the second example (n = 2) this trend is reversed. 
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We now proceed to study the normal modes of radial pulsations of the equilibrium 
configurations. To this end we explicitly compute the solutions to the boundary value 
problem consisting of ( 1271) with the boundary conditions u(0) = and u(R) = const. 
This task is technically nontrivial because with the EoS (1311) the point r = R is a 
regular singular point of (1271) . as can be shown by the following consideration. The 
power expansion of the anisotropic TOV (133]) reveals that as r — > R, 



Po^l ^l ■ ^ H - (35) 



2k(n + - fi s ) J V R 

where // s = 2M/R is the surface compactness. With the above result the asymptotic 
behaviour of the coefficients of (126 al) - fl26 c|) as r — > R follows as 

(1 - /i s )/i s ( r\ ( \\ (1 -/i s )/i s (8-7/x s )/i s 

^^^^l 1 -^' Cl ^"^ + nJ^^' C ° 4iP • (36) 

Due to the vanishing of C*2, the point r = i? is a regular singular point of ( 1271) . which 
allows for unbounded solutions. However, regularity of the solutions can be enforced 
by integrating (|27|) as an initial value problem starting from r = R with the initial 
condition u(R) = const, and 

u(R) = ~ u(R) = -u{R). 37 

d 2R{1 + l/n)(l - // s ) 

(We note that due to simplicity of the EoS ( 13T|) the above boundary condition coincides 
with equation (7c) of (13] that was derived assuming an asymptotically polytropic EoS 
describing the fluid with isotropic pressures.) Due to the robustness of the numerical 
method [29], finding the eigenvalue uf and the eigenfunction «j(r) satisfying the above 
boundary condition and also w(0) = with i = 0,1,... nodes in (0, R) for any 
equilibrium configuration is now a simple task. We are most interested in i = since, as 
discussed in the preceding section, if > the equilibrium configuration is considered 
stable, while if < it is considered unstable. 

For all solutions shown in figure [TJ we computed the eigenvalue and observed 
that as the central energy density increases (n, k and a fixed) it changes the sign from 
positive to negative exactly where, within the remarkable numerical precision of the 
procedures, the M(R) curves exhibit their maxima. In other words, the two methods we 
have applied to assess the stability properties of the equilibrium configurations obtained 
from the quasi-local EoS gave compatible results. 

In figure [2] we show a highly compact configuration obtained with n = 1, k = 
100 km 2 and a = 4, and central energy density for which Uq = (more details can 
be found in the figure caption.) This configuration is therefore considered marginally 
stable. An interesting feature that we can observe is the energy density which is not 
maximal in the centre of the configuration, but is outwardly increasing over the central 
region. 
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Figure 2. Critical (marginally un/stable) configuration with parameters k = 100 km 2 , 
?i = 1, a = 4, p (0) = p c ^ 1-87 x 10 15 gcm- 3 , M ~ 2.14M©, i? ~ 9.14km, 2M/i? ~ 
0.69. Energy density (thick), radial (long-dashed) and transverse pressures (short- 
dashed), and the first three eigenfunctions (thin lines), eigenfrequencies corresponding 
to the first three eigenfunctions are Uq ~ 0, v\ ~ 4.75 kHz, ^2 — 6.89 kHz, respectively. 



5. Summary and outlook 

We have shown that with the quasi-local variables one can construct equations of state 
(EoS) to describe the anisotropic fluid in spherical symmetry. The quasi-local EoS may 
be used to generate static anisotropic spheres over a wide range of radii and masses, by 
varying the central energy density as the parameter. We have also studied the stability 
of the configurations with respect to radial perturbations. The analysis of normal 
modes of radial pulsations and the M(R) method gave identical results. The numerical 
analysis of several examples revealed that when the transverse pressure is greater than 
the radial pressure, the maximal surface compactness (2M/R) of the spheres that can 
be obtained before the onset of instability is higher than that of isotropic spheres. 
Another interesting feature is the outwardly increasing energy density in the cores of 
some stable configurations, despite of the fact that the radial pressure and the energy 
density are related by the polytropic EoS (figure [2]). (Some may recall Weinberg's 
remark on isotropic spheres that "it is difficult to imagine that a fluid sphere with a 
larger density near the surface then near the centre would be stable" [30] ; we saw here 
that a similar structure can be supported by anisotropic pressures, whatsmore, it can 
be stable.) Therefore, the concept of the quasi-local EoS has been shown to be useful 
in constructing interesting models of astrophysically plausible objects such as highly 
compact stars, especially since the M(R) method can be used to assess their stability 
properties. 
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